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We exhibit measurements for optimal state estimation which have a finite number of outcomes. 
This is achieved by a connection between finite optimal measurements and Gauss quadratures. The 
example we consider to illustrate this connection is that of state estimation on N qubits, all in a 
same pure state. Extensions to state estimation of mixed states are also discussed. 
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Two central questions in quantum information theory are: "How much information can be encoded in a finite 
quantum ensemble?" and "What measurements should one use to optimally retrieve this information?". There are 



^vq ■ many ways to address these questions. One of the simplest is that proposed in 0, in the form of a game between a 



referee, Alice, and a player, Bob. Alice prepares N qubits in an identical state |0), unknown to Bob. This state is 

drawn from a uniform probability distribution on the set of all pure qubit states, the so-called Bloch sphere. Then 

Bob performs some measurement on 10)®^, to guess the prepared state |fi). If the state guessed by Bob is |flfc), he 

achieves the score | (Olfi/.) | 2 . Bob aims at achieving the best possible average score 
00 . 

Csl ■ s w = J2 fdnp N (n k \n) \(n\n k }\ 2 , (i) 

> : k J 
r» . 

£C) • where J dft = 1 and where Pjv(fife|f2) denotes the conditional probability of Bob to get a measurement outcome |Qfc) 
CN | when Alice has prepared the state Ifi)®^. 

It is proven in fij, that when N identical copies of the qubits whose state we want to estimate are available, the 
optimal value that Eat can take is 

■ 

O ' AT + 1 

Z^UN) = — . (2) 

^ 1 V ' N + 2 W 



Also, one readily checks that the following positive-operator-valued measure (povm) achieves the optimal score: 



§■ {N + i)\n)(n\® N dn. (3) 

I n t ne sequel, we will call this povm the coherent povm, because it is made of projectors onto the so-called SU(2) 
coherent states Q- While this povm requires an infinite number of outcomes, it is possible to build povm's achieving 
the optimal score S 1 ^* and whose number of outcomes is finite. This problem has been addressed by several authors. 
In an algorithm to derive finite povm's achieving the score @ was proposed. Later, in Q|, the issue of finding 
H | minimal optimal finite povm's was tackled and a solution was provided for up to 7 qubits. 

The present note aims at bringing a new perspective on this problem. We show that there is a natural connection 
between this problem and the design of a Gauss quadrature on a sphere (see also 0). Let us first introduce some 
notations and reformulate the game between Alice and Bob in more precise terms. In the following, {|0), |1}) will 
denote an orthonormal basis of the Hilbert space of a qubit. Pure qubit states will be parametrised as 

|O) = cos^|0)+e^sin^|l). (4) 

where < 4> < 2ir and < < ir. CI = (6, cf>) denotes the direction on the Bloch sphere associated with the qubit 
state. 

The smallest Hilbert space describing the states |J7) (2l7V is the symmetric subspace of the Hilbert space of N qubits, 
3{f N (that is the subspace of maximal total spin N/2). As any measurement on !K® N , the measurement used by Bob 
is a povm describable by a set of operators {0^} satisfying 

fe > 0, (5) 
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J2°k=lN, (6) 



k 



where Ijv denotes the identity on . 

Now observe that any set of n operators Ofe = c k \^l k ) l^l k \® N such that 



J2ckPk){n k \® N = iN- (7) 
fc=i 

is an optimal povm. Indeed this povm is such that P/v(^fc|^) = Ck\(Qk\ty\ 2N , it achieves the score 

z N ({o k }) = ^c fe y dn|<n fe |n>| 2 ^ +1 >. 



/c=i 

As a consequence of Shur's lemma, we have 

dn\n)(n\® N+1 = ^±L. 
i /\ i N + 2 

Thus the score achieved by the povm J7|) reads 
1 . - N 4- 1 

^({o fe}) = — E^ = FT^ = E - (8) 
fc=i 

The issue of finding optimal finite povm's can actually be re-stated in a very simple form, that allows for a straight- 
forward derivation with n ~ (N + l)\ (N + 1 ) /2] for arbitrary values of N , as well as other more " economical" optimal 
povm's for some large N values. Our construction is based on the concept of Gauss quadrature which we briefly 
describe now. 

A Gauss quadrature is a rule designed to approach an integral over some domain D by a sum and defined by a 
mesh of points x k € 2) and associated weights Aj, , 

f « Y>/(z fe ). (9) 

If the parameters {cc/-, \ k } are chosen so that the quadrature is exact for a set of functions fj(x): 
[ dx.f j (x)=J2x k f j (x k ), (10) 



it will of course be exact for any linear combination of the form f(x) = ■ a.jfj(x). 

In order to associate a Gauss quadrature with our problem, let us observe that Eq.Q can be expressed as 

n 

j^ckWkm 2 = i vn. (ii) 

fe=i 

Expanding (f2fc|f2) in terms of spherical harmonics YJ m (f2), this condition can be re-expressed as the following 
system of equations 4] : 



J2c k = N+l (12a) 
k=l 

n 

Y^ckYr^ik) = o, (12b) 
fe=i 

where I = 1, . . . , N; m = —I, . . . , I. Using the following properties of spherical harmonics: 

Y°(n) = l/yfe, (13) 
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J dSlYT(Q)Y^'(Q) = S w 6 mm ,, (14) 
we see that the conditions (|TT|) can be rewritten as 

~ n 

dciYr(ii) = j2^Yr(^ k ), (15) 
k=i 

where X k = 4irc k /(N +1). Thus Ea. l|15|l means that finding an optimal finite povm amounts to finding an exact 
Gauss quadrature on the sphere for the spherical harmonics up to order I = N. 

We now show that some quadratures can be found very simply. Spherical harmonics can be written as: 

Y, m (fl) = C/"P Z H (cos 9)e im ' 1 ' (16) 
where C™ are constants (whose precise expression are irrelevant here), 

p r (u) = (l-u^—P l (u) 7 (17) 
are the associated Legendre functions and 

w = ^(. 2 -iy as) 

are the Legendre polynomials. We see that the spherical harmonics are either a product of e lm ^ and a polynomial of 
degree / in cos# (when to is even), either a product of e 1 " 1 ^, sin# and a polynomial of degree I — 1 in cos 9 (when to 
is odd). 

Thus, expressing Eq. l|15|l in spherical coordinates, we must find a Gauss quadrature that is exact for: 

d(cos9) / d<pcos n 6e im 'l' (19) 
1 Jo 

V to even (0 < in < N) and n (0 < n < N), and exact for 

1 p2tt 

d cos 6 / d</)sin6)cos n 6<e im0 (20) 
-l Jo 

V to odd (0 < to < N) and n (0 < n < N - 1). 

A simple Gauss quadrature is obtained by separating integrations over the variables 6 and 4> an d defining an 
independent Gauss quadrature for each. Let us first consider the integral in ip. Since to is integer, we have: 

d<f> e im * = 6 m o (21) 

Since to can take integer values up to AT, it is easy to see that the minimal Gauss quadrature that will be exact 
for these N + 1 integrals will require n\ = N + 1 equidistant points, for instance the (iV + l) th roots of unity 
4>i = ^l^TT (0 — •? — -^0' eacn w ith a same weight t|tT- 

Let us now consider the integral in 9. Ea. (|21|l shows that integrals for to ^ will cancel out if the quadrature for 
</) is exact. So, the quadrature for 9 needs only be exact for to — 0, that is for the integrals 

d(cos 9) cos" 9=1 dxx n V < n < TV. (22) 
-l J-i 

The following theorem [6j provides us with such quadratures: 

Theorem 1 Let [a, b] C R, and let {p n } denote a complete set of orthogonal polynomials on L 2 ([a, b]). If X\ < .. . < 
x n denote the zeros of {p n (x)} ; there exist real numbers Ai, . . . , X n such that 

/b n 
da(x)p(x) =^2\ k p(x k ), (23) 
fe=i 

whenever p(x) is an arbitrary polynomial of degree at most 2n — 1. Moreover, the distribution da(x) and the integer 
n uniquely determine these numbers \ k - 
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Using this theorem allows straightforwardly to build a Gauss quadrature on [—1, 1] using ri2 = \(N + l)/2] points, 
that is exact for any polynomial up to degree N, such that we finally obtain a povm with 



ni n 2 = (N + 1) 



N+l 



(24) 



elements. As the #-part of the quadrature is based on Legendre polynomials, this separated variables Gauss quadrature 
is generally known as Gauss-Legendre quadrature. 

The derivation of optimal povm's from the Gauss-Legendre quadrature is generic but we have no guarantee that 
it provides us with a minimal optimal povm. Actually, it is not the case, as we already know that Latorre et al 
have been able to construct optimal finite povm's with less elements than ours, at least for N < 7 One could 
have expected that the Gauss-Legendre quadrature we have presented above is not minimal. Hcuristically, finding a 
minimal Gauss quadrature is equivalent to finding n points on a sphere so as to optimally cover it. Intuitively, one 
expects that minimality will in general not be achieved by considering the 2 variables 9 and <$> separately as we did. 
In a series of articles 0H,0, Lebedev deals with the issue of defining "good" quadratures and provided a quadrature 
rule on the sphere for all odd values up to N = 29, as well as for higher values in the form iV = 6a + 5 up to iV = 131 
(or a = 21). This last sequence requires a number of points 

(A + l) 2 , , 

n = 3 ' + 2, (25) 

improving on Gauss-Legendre (separated variables) quadrature and on the conjectured minimal number of points 
given by Eq. (24) of g . 



FIG. 1: Lebedev quadrature of order N = 59. Distribution of 1800 points for Lebedev quadrature of order N = 59 (taken from 
0). 

Fig- El shows one of the quadratures proposed by Lebedev in 0. As may be seen, quadratures may be "polarised", 
in the sense that not all points are equivalent by symmetry, moreover not all points are associated with a same weight. 

A uniform distribution on the sphere is actually achievable for some number of points n, and is called a spherical 
iV-design if the Gauss quadrature it defines is exact for all spherical harmonics up to order N: 

/A " 
dfl Y{ n {n) = —J2 Y " l ( n k)J = 1 ... iV, ra =—/... Z. (26) 
71 k=l 

Spherical iV-designs are known up to N = 13 10]. An example is displayed in Fig. [3 

A comparison between the various quadratures we have discussed is shown in the four first columns of Table Q] 
The previous results can be straightforwardly applied to mixed states. Suppose that we are given N copies of a 
qubit mixed state p(rj, n) = ^(I + r\ n • a), where the orientation n is drawn from a uniform distribution over the 
2-sphere, and where the shrinking factor 77 is drawn an arbitrary probability distribution f(rj) over the interval [0 : 1]. 
The problem of estimating at best the state p(rj,n) has been addressed in A key result of is a prescription 
of optimal state estimation strategies for mixed states from the knowledge of optimal state estimation procedures 
for pure states. Therefore, the povm's we have derived in the previous sections can immediately be used to obtain 
finite optimal povm's for mixed states. A second result of [Tl). which is relevant here is a formula for the minimum 
number of povm elements for optimal mixed state estimation, n^^. Interestingly, this formula is independent of the 
distribution f(rj). It reads: 



(AT) 
mill 



N/2 
s=so 



(27) 



FIG. 2: Spherical design with n = 60 points defining a Gauss quadrature on the sphere of order N = 10 (lo||. 



where sq — if N is even, and sq = 1/2 if N is odd, n(0) — 1 and n(s) is the minimal number of povm elements for 
pure state estimation on 2s qubits. From the pure state povm's based on the Gauss-Legendre quadrature, we find 
the (loose) upper bound < Ef= 7 s 2 (2s + l)\(2a + l)/2], from which we find 

(N) N 3 5 9 17 , x 

<il < -^+8 N + U N + 1 ( Neven )' ( 28 ) 

n^<^ + ^ + y 2 N + \ (Nodd). (29) 

Furthermore, for odd values of N up to 29, we may improve on this bound by using the Lebedev quadrature, as shown 
in the last column of Table [I] 

Finally, let us show that this approach also allows to deal with extensions of this problem, such as the estimation 
of general states in the tensor product of the N Hilbert spaces 0], instead of states restricted to the symmetric 
subspace J{® N . More precisely, suppose we want to perform state estimation on a set of states 

{U^ N (gM)^\^g)),geSU(2)}, 

where now \<j>) is not restricted to belong to the symmetric subspace 3~C® N . Rather, we have 

M (30) 

i 

where the index i runs over the irreducible representations of U® N . Accordingly, U® N (g) = Yli^i(9) anc ^ I^CsO) = 
ai\4>i(g)) , where \(f>i(g)} = Ui{g)\4>i). Suppose now that we have an optimal covariant povm for this task: 

{U{ 9 r N uHgf N = 0(g),ge SU(2)}, (31) 

The score is given by 

S= [ dg x ( rfff 2 |(3i|32)| 2 (0(5i)|O(52)|0(5i)) 



= JdgiJ dg2\(gig2 1 \e)\ 2 (cP(gig2 1 )\0\<P(gig2 1 )). 

With the change of variable g = gig^ 1 and the fact that J dgi = 1, we get 

5 = | dg\(g\e)\ 2 (Hg)m(g)). 

Again, any set of positive operators {ckO(gk)} that satisfies J2k c k®(dk) = X)i I* i s an optimal povm (I; denotes 
the identity over the irreducible subspace Since %i = {\(f>i(g)),g S SU(2)}, this condition is equivalent to the 

following 

£cfc<&Cff)|Q(fflOl0i(ff)) = 1 Vfl.i. 
k 
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Let us diagonalize the povm elements as = l0 (/3) )(0 (/3) l, and define the 577(2) elements # (/3) such that |0 (/3) ) 
decomposes into the subspaces JCi as follows 

i 

We may now write the povm condition as 

k 

which will be satisfied as soon as 

J2^\(Mg)\M9k))\ 2 = 1 Vff.t. (32) 

We therefore have a condition equivalent to Eq. (11) for each subspace Since each "Ki corresponds to a space of 
total spin Si < N/2, we need to find a Gauss quadrature on the sphere that is exact up to order 25*1 + 1 for subspace 
JCj. It therefore suffices to consider the largest subspace, that is the symmetric subspace 3~Cf N , corresponding to a 
maximal total spin Si — N/2, which requires an exact Gauss quadrature up to order N + 1, and taking the same 
quadrature for the smaller subspaces will solve the problem. Our solution for the symmetric subspace therefore 
immediately yields a solution for the whole space, using a povm with the same number of elements. 

In conclusion, we have established a connection between the issue of finding optimal povm's and the problem of 
finding Gauss quadratures on a sphere. Following this idea, we easily derive a solution for arbitrary large N, that while 
not minimal, already improves on the best known solution Q. Moreover, with the help of mathematical literature on 
Lebedev quadrature [3, la, |9j or spherical designs [10( , we are able to give enhanced solutions for some large values 
up to N — 131. Although the connection doesn't solve the problem of finding minimal optimal povm's, at least it 
helps to understand the difficulties involved since they turn to be essentially the same as that for finding good Gauss 
quadratures. Apart from its theoretical interest, this connection could be used to design optimized measurements for 
experimental state estimation. Even though the exact implementation of these povms for large number of copies N 
would probably prove to be hard to realize experimentally, they could be used as a starting point to derive feasible 
approached (non- minimal or non-optimal) povms (see also the related works by Hayashi [13j, and Gill and Massar 
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N 


Latorre a/ 


Legendre 


Lebedev 


spherical design 


Legendre (mixed) 


Lebedev (mixed) 


1 


2 


2 




2 


2 


2 


2 


4 


6 




4 


7 




3 


6 


8 


6 


6 


10 


8 


4 


10 


15 




12 


22 




5 


12 


18 


14 


12 


28 


22 


6 


18 


28 




24 


50 




7 


99 


■ >— 


9fi 


94 


fin 




8 




45 




36 


95 




9 




50 


38 


48 


110 


86 


10 




66 




60 


161 




11 




72 


50 


70 


182 


136 


12 




91 




84 


252 




13 




98 


74 


94 


280 


210 


14 




120 






372 




15 




128 


86 




408 


296 



TABLE I: Comparison of the number of elements n of optimal povm's built using different methods as a function of the number 
of copies N of the estimated state. The first column corresponds to the minimal solution proposed by Latorre et al in A\ up to 
N — 7. In the second column, we see that the separated variables Gauss-Legendre quadrature provides an optimal povm for 
arbitrary large N, although the number of elements N ~ (N + 1) 2 /2 is not minimal. The third column corresponds to Lebedev 
quadrature that reduces the number of elements to N ~ (N + l) 2 /3 but only exists for some odd values up to N = 131. The 
spherical designs, that define povm's with equally weighted elements and are known up to iV = 13, are given in the fourth 
column. The last two columns contains a bound on the minimal number of povm elements necessary for mixed state estimation 
using Gauss-Legendre quadratures and Lebedev quadratures. 



